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ABSTRACT 

Context. Various morphologies such as multi-concentric shells and spiral-like patterns have been observed around Proto-Planetary 
Nebulae and AGB stars. It is widely argued that the regular density patterns are produced by binary systems. 

Aims. The goal is to build up a simple dynamical model for the out-flowing circumstellar envelope around AGB stars in a wide binary 
system to explore the parameter dependence of the geometrical characteristics of column density patterns. 

Methods. For an AGB star in a wide binary system, the orbital motion of the star can be approximated as a series of pistons that 
simultaneously push out dust and gas along radial directions, but work in different oscillation phases. The piston model can fast 
produce column density patterns with high enough spatial resolutions for parameter dependence exploration. 

Results. The formation of 3-D quasi-spherical density structure is induced by orbital motion of the AGB star. The column density 
pattern only depends on two parameters: eccentricity of the orbit e and the terminal outflow velocity to mean orbital velocity ratio 
y. When viewed perpendicular to the orbital plane, spiral, broken spiral, and incomplete concentric shell patterns can be seen, while 
when viewed along the orbital plane, alternative concentric half-shell, egg-shell, and half-shell half-gap patterns will develop. Non- 
zero eccentricity causes asymmetry, while larger y makes a weaker pattern and helps bring out asymmetry. A spiral pattern may 
becomes broken when e > 0.4. The spiral center is always less than 12% of spiral pitch away from the orbit center. One should have 
more chances (~ 80%) seeing spiral-like patterns than seeing concentric shells (~ 20%) in the circumstellar envelope of wide binary 
AGB stars. 

Key words. Stars: AGB and post-AGB - Stars: circumstellar matter - (Stars:) binaries: general - shock waves 



1. Introduction 

In recent years, various kinds of morphology in circumstel- 
, £^ ' lar envelopes of AGB stars have been re vealed by high spa- 
k> \ tial r esolution observations. For example, iMauron & HugginsI 
•jj ■ (I2000h showed a series of shells in the envelope o f the carbon 
d ' rich A GB star IRC +10216 in an optical image. Bali ck et"ai] 
d2001l) showed clear multi-shell structures in the relic AGB en- 
velope arou n d the planetary nebula NGC 6543, the Cat's Eye. 
iKwok et alJ §998) found multi-arcs near the two lobes of a 
bipolar proto-planetary nebula IRAS 17150-3224. High veloc- 
ity bipolar jets had also been discovered in the AGB circumstel- 
lar envelope by mea ns of interferometry of maser emission, e.g., 
the observations by (llmai et alj2002ll2005h towards W 43A. The 
most mar velous pattern found most r ecently is the spiral pattern 
found by Mau ron & HugginsI d2006l) in the envelope of AFGL 
3068, a C-rich AGB star. They fitted an Archimedes' spiral to the 
observed spiral pattern, which demonstrated the amazing regu- 
larity of the pattern. 

The mechanis m that produces the various m orphologies is 
still controversial. Steffen & Schonberner (2000) proposed peri- 
odic change of mass-loss rate or wind speed, but the change of 
mass-loss ra t e by th ermal pulse has too long a period (~ 10 5 yrs). 
Bali ck etail (|2P01) attributed the multi-shells around the Cat's 
Eye to instability of physical and chemical processes during 



dust formation and growth, as suggested by Sim is et al ] (l2001h . 
who tried to ex plain the concentric shells a round IRC +10216 in 
the same way. IMauron & HugginsI d2006l) preferred to explain 
their spiral pattern in AFGL 3068 by numerical results from 
Mastrode mos & Morrisl dl999l) . (MM99) hereafter, in which an 
AGB star in a wide binary produces a spiral-like density pat- 
tern by the reflex motion of the star (their model 4). MM99 
performed smooth particle hydrodynamic (SPH) simulations to 
explore the effects of binarity on AGB wind for circular orbit 
binaries. In their wide binary example, they showed spiral and 
concentric shell density patterns formed in the circumstellar en- 
velope. But to understand the observed patterns better, it's help- 
ful to know, if the orbit is not circular, what will be different 
in the patterns. However, due to difficulties in SPH simulation, 
e.g., time consuming, not high enough spatial resolution, vary- 
ing of the particle's smoothing length, etc., a simpler simulation 
is needed for extensive parameter exploration tasks. 
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ISokerld 19941) discussed the influence of a wide binary on the 
structures of planetary nebulae by taking the AGB mass-loss as 
a sequence of mass loss pulses. With far distance piston approxi- 
mation and sticky particle model in my work, the influence of bi- 
narity will be explored in more detail. Section|2]gives the numer- 
ical piston model to simulate the formation of density pattern in 
the circumstellar envelope of a wide binary AGB star. Numerical 
results are present in Sect. [3] A brief summary is given in Sect. [4] 
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2. Piston model and column density pattern 

2.1. Piston model 

To simplify the problem, several assumptions about an AGB star 
are made: mass-loss rate does not change with time; outflow ve- 
locity is constant and equal to terminal velocity V e ; gas pres- 
sure can be neglected; and dust gas interaction can be neglected. 
Although the mass-loss rate of AGB stars may change with time 
due to thermal pulse, it chan ges on a very long timescal e (e.g., 
10 5 years, as suggested by Steffen & Schonberner 2000) that is 
much longer than the orbital period. Mass-loss rate and outflow 
velocity m ay change due to the interplay of dust formation and 
dynamics ( Si mis et"aT]|2001l) . but it is not clear whether such a 
phenomenon really occurs. The assumption of constant outflow 
velocity is equivalent to assuming that the gravity force and ra- 
diation pressure on dust grains and gas are balanced everywhere. 
This is supported by hydrody namic calculations (e.g., the model 
results for IRC +10216 bv lSkinner et alJI 19991) . unless one ex- 
amines the region very close the star. Gas pressure can broaden 
any density structure by means of thermal dissipation, but its ef- 
fect is marginal in the outer part of the supersonic outflow. And 
gas pressure does not alter the geometry of density structures, 
which is the main topic of this work. Gas and dust are not distin- 
guished and the material is always called gas. These assumptions 
sound plausible for this work because hydrodynamic effects on 
the pattern geometry become unimportant in wide binary sys- 
tems. 

Under these assumptions, the motion of each chunk of 
ejected gas can be treated as ballistic motion. One can imag- 
ine that many gas parcels are ejected from the star surface with 
constant velocity V e relative to the star. The absolute velocity of 
each gas parcel in the center of mass (CoM) frame of the binary 
is equal to the sum of V e and star motion velocity. 

This work only considers density patterns formed far away 
from the star in the circumstellar envelope. In this case, the star 
motion and mass loss can be approximated as a piston that is 
ejecting gas parcels and moving back and forth along radial di- 
rection to modulate the velocity distribution in the ejected gas. 
The ejected gas will be compressed by the velocity gradient 
to form a density enhancement. Shock may happen when gas 
parcels collide with each other due to different velocities. One 
such piston works in one direction. The whole three-dimensional 
AGB mass-loss envelope can be imagined as a collection of 
pistons pointing in different directions and running in different 
phases. The different piston phases in different directions con- 
catenate the density enhancements or shocked regions to form 
a three-dimensional quasi-spherical density structure, like the 
structure produced by the SPH simulations by MM99 for cir- 
cular orbit in their model 4. This is called a far distance piston 
approximation. 

The orbital motion of the AGB star is assumed to be 
Keplerian. Projection of the star motion along any given direc- 
tion (9,(p) (see definition of quantities in Fig. [TJ gives the posi- 
tion of the piston along that direction as 



l(t) = a ? (cos £p cos 4> + Vl - e 2 sin E v sin <p} sin 6. 



(1) 



Here e and a p are the eccentricity and semi-major axis of the star 
orbit, respectively, while E p is the eccentric anomaly of the star. 
Then the speed of the piston is 



Fig. 1. Schematic illustra- 
tion of the orbital motion 
of AGB star (point P), with 
definitions of Cartesian co- 
ordinates OXYZ, eccentric 
anomaly E p , and position 
angle (#,</>) of the piston 
(the heavy arrow). Point O 
is the center of the orbit 
while G is the center of 
mass. The dotted circle is 
the circumscribed circle of 
the orbit. 



where a> — 2n/T is the mean angular velocity of the star motion 
with orbital period T. The time dependence of piston motion is 
reflected by 




E p - e sin E p 



u)t. 



(3) 



In the numerical simulation, the An solid angle is divided 
into 10000 equally sized patches and a piston is computed along 
the center direction of each solid angle element. Each piston is 
assumed to eject 100 gas parcels (with equal mass) within one 
orbital period. Gas parcels ejected at different times set out from 
different positions l(t) with different velocities V e + Vi(f). After 
leaving the star, each gas parcel moves outward with constant 
speed, until it collides with neighboring parcels. The collision 
of gas parcels causes shock and the collided gas parcels are ar- 
bitrarily assumed to merge with each other (actually move to- 
gether side by side with very small fixed distances in the simu- 
lation, sticky particles). The velocity of the merged gas parcels 
is determined by conservation of momentum. The spatial distri- 
bution of all gas parcels represents the 3D density distribution of 
out-flowing material in the AGB envelope. Projection of the 3D 
distribution of gas parcels onto any specified sky plane produces 
a column density pattern. 

2.2. Normalization of the problem 

Several physical parameters are involved in the piston model: 
orbital parameters such as semi-major axis a p , eccentricity <?, and 
period T, and the steady outflow velocity V e . However, not all of 
them are independent of the problem. To reduce the number of 
parameters, the above formulae need be normalized by dividing 
all length quantities by the circumference of the circumscribed 
circle of the orbit 2na p , dividing all time quantities by period T, 
and dividing all velocity quantities by the mean orbital velocity 
of the star wa p . Then, the piston motion formulae Eqs. (Q]i and 
(O become 

l(t ) = — (cos Ep cos <p + Vl - e 2 sin E v sin (p) sin 6 (4) 



V 1 - e 2 cos E D sin d> - sin E B cos i 
V/(f) = ^ - sin 



1 - e cos E„ 



and 



E p - e sin E p = Int. 



(5) 



(6) 



dl(t) Vl - e 2 cos E D sin <b - sin E„ cos <f> 

V,(t) = = coa p — ? — — P - sine, (2) 

df 1 - e cos E-o 



In the numerical simulation, each gas parcel is ejected at nor- 
malized position 1(f) with the constant normalized speed y+ V/(f) 
(with y = Ve/wa p being a velocity ratio). After the normaliza- 
tion, one has only two free model parameters: e and y. All other 
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parameters such as a v , T, and V e are combined in the parameter 
J- 

2.3. Formulation of single piston density pattern 

Thanks to the simplicity of the piston model, it is possible to 
directly formulate the density pattern geometry (position of den- 
sity peaks or shocked regions) in the outflowing gas. In some 
cases, solving the formulas may be easier than analyzing numer- 
ical column density maps in the discussion of certain aspects of 
the pattern geometry. But due to the limitation of space, the de- 
duction of normalized pattern formulas is put in AppendixlAlthat 
is available only in electronic version. 

3. Results 

3.1. Numerical column density patterns 

Column density pattern geometry has been explored in a large 
range of parameter space: e = 0.0 - 0.8, y = 10 - 90, and dif- 
ferent observers' direction angles (£, tfr). Representative types of 
morphology found in the computed column density maps are 
shown in Fig. [2] In the first row patterns are viewed perpendic- 
ular to the orbital plane, while in the second row patterns are 
viewed along the orbital plane. Maps A and D show that, for the 
circular orbit case, the column density pattern is a good spiral- 
like pattern when viewed perpendicular to the orbital plane and 
alternative concentric half shell pattern when viewed along the 
orbital plane, which closely resembles the SPH simulation re- 
sults of MM99 (their Model 4). Comparing map A with B tells 
us that, when eccentricity is non-zero, the spiral pattern may be- 
come broken in its inner part (the broken spiral pattern in map 
B). The extent of the spiral break decreases when the shock pat- 
tern moves outwards, and eventually continuous spiral forms ap- 
pear in the outer part of the pattern. If both e and y are large, 
as the case of map C, even if observed perpendicular to the or- 
bit plane, the column density pattern still looks like a series of 
concentric shells, but the shells are not only much fainter, but 
also incomplete in the lower-right quarter {incomplete concen- 
tric shell pattern). 

A common character of the patterns in the second row of 
Fig.|2]is concentric shells with the lack of density pattern along 
the vertical directions, which are the directions of the orbit pole. 
Comparing maps D and E tells us that non-zero orbital eccentric- 
ity may cause the variation of relative radii of the concentric half- 
shells and make strong non-shocked patterns around shocked 
patterns. As a result, in map E, the radii of the left and right 
series of half-shells are nearly the same, which imitate nested 
egg shells. In the right half of map E, the innermost shells look 
thicker. They are actually double shells with non-shocked gas 
filled in between. Moving outwards, the double shells eventu- 
ally merge with each other to form sharp single-layer shocked 
shells. One should also have noticed the left-right asymmetry 
both in density contrast and in pattern shape in map E, which is 
the effect of non-zero eccentricity, too. Comparing map F with 
maps D and E gives us a surprise: the larger e and y in map F 
cause the missing of shells on the right side, which results in a 
half-shell half-gap pattern. 

Exploration of more maps shows that left-right asymmetry 
(as in the broken spiral, incomplete concentric shell, egg shell, 
half-shell half-gap patterns) only appears when e > and y 
is large, otherwise, when y is small, even if e > and the 
asymmetry exists, the asymmetry may be buried in the invisible 
central region of the maps. When the asymmetry is prominent, 




Fig. 5. The border line di- 
vides the y — e parameter 
space into two parts: the 
density pattern in the or- 
bital plane of the model is 
a continuous spiral pattern 
in the lowerleft part and a 



broken spiral or incomplete 
10 20 30 40 so 60 70 8o 90 100 concentric shell pattern in 
' the upperright part. 



its dependence on the observer's azimuth angle is also found 
when observing along the orbital plane: the pattern looks more 
left-right symmetrical when being viewed along the major orbit 
axis, while it appears more left-right asymmetrical when being 
viewed along the minor orbit axis. Several maps that illustrate 
this effect are present in Fig.|4]in the electronic version only. 

3.2. Several interesting properties of the column density 
patterns. 

One of the interesting problems is how the transition from spiral- 
like pattern to concentric half shell pattern happens when the 
observer's inclination angle changes from 0° to 90°. In Fig. [4] 
more column density patterns are shown for different observa- 
tion inclination angles f = 0°,68 o ,79°, and 90° for the circular 
orbit case. One can see that only when the inclination angle is 
f > 79°, does the column density pattern begin to prominently 
deviate from the spiral-like pattern. That is to say, in a large 
range of inclination angles, from 0° to 79°, one can see a spiral- 
like pattern. If one sees a concentric shell pattern, the inclination 
angle of observation must be within the 90 ± 10 degree range, 
that is, the star orbit must be almost edge on. If the distribution 
of the star orbit orientation is completely random, one has an 
approximate probability of 1 - cos 79° =81% seeing the spiral- 
like pattern and only a probability of 19% seeing the concentric 
shell pattern in the circumstellar envelope of such a wide binary 
AGB star. MM99 intuitively speculated that, when viewed with 
finite resolution at great distances from the star and from some 
intermediate latitude, the multi-shells may appear to form full 
circles. This is contradicted with the conclusion of this work and 
is incorrect. Actually the relative shift of shells with inclination 
angle is true but the shift of shells forms a spiral pattern, instead 
of full rings. This can be easily understood because one sees spi- 
rals when observing perpendicular to the orbit plane and because 
the left and right shells have different radii. 

Another question is when the spiral pattern becomes a bro- 
ken spiral. The answer lies in the single piston model along the 
pericenter direction of the orbit because the broken part of a spi- 
ral always happens on the pericenter side. Some trial solutions 
of single piston formulas dA.41 IA.51 IA.6I ) demonstrate that the 
formation of the broken spirals is due to a double solution of E p 
within one period. These formulas are numerically solved for 
all possible combinations of (<?, y). Those pairs of (<?, y) beyond 
which a double solution of E p begins to appear are shown in 
Fig.|5]as a border line that divides the parameter space into two 
halves: lower left half for the continuous spiral; upper right half 
for the broken spiral and the incomplete concentric shell pattern. 
Because e does not change too much around 0.4 along the border 
line, one may say that a broken spiral pattern may appear when 
eccentricity e > 0.4. 
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Fig. 2. Column density pat- 
terns calculated by the 
piston model for repre- 
sentative cases. Observer's 
azimuth angle is always 
6 = -90°, which means 
the pericenter of star or- 
bit is pointed to the right. 
Inclination angle £ and e, y 
are shown on top of each 
map. All maps are in lin- 
ear stretch in the same gray 
scale range. X and Y axes 
are in normalized length 
unit (i.e., times of 2na p ). 
Pattern types shown at the 
right bottom of each map 
are: sp - spiral; bsp - bro- 
ken spiral; ics - incomplete 
■ circular shells; acs - alter- 
(hshg) na (j ve concentric shells; es 
- egg shells; and hshg - 
half-shell half-gap. 



C=79° 




■60 -40 -20 2D 40 60 
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-B0 -40 -20 20 40 
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Fig. 4. Maps show inclina- 
tion angle effects. Azimuth 
angle is t// = -90°, e = 0, 
and y = 10 for all maps. 
One can see the spiral pat- 
tern in the inclination an- 
gle range of 0° - 79° and 
the concentric shells in the 
79° - 90° range. 



The third question is where the spiral center is. Archimedes' 
spiral curve is fitted to the mature part of shocked patterns com- 
puted from the piston model in a wide range of parameter space 
(e = - 0.8,y = 2 - 90). The pattern fitting demonstrates that 
the shocked patterns are always good Archimedes' spiral curves. 
The center of the spiral is not at the orbit center, but the shift 
of the spiral center from the orbit center shows some degree 
of dependence on both parameters e and y. Generally speak- 
ing, higher eccentricity and smaller y causes larger spiral center 
shifts. However, in most cases, the spiral center shift is always 
less than 2% of a spiral pitch in the major orbit axis direction and 
always less than 12% of a spiral pitch in the minor orbit axis di- 
rection, which means the spiral center should be just on the star 
position when the observation resolution is not extremely high. 
In circular orbit cases, the spiral center is exactly at the orbit 
center. 



The properties of the various column density patterns may be 
critical in distinguishing the binary mechanism from the mass- 
loss rate modulation mechanism. The patterns also serve as a 
magnifier to show the binary motion in the center of the circum- 
stellar envelope. On the other hand, the formation of the density 
pattern, especially the shocked density pattern, greatly enhances 
the visibility of the circumstellar material and enables us to study 
the mass-loss history of AGB stars in binary systems. 



4. Summary 

A piston model is used to simulate the quasi-spherical density 
patterns formed in the supersonic outflow of wide binary AGB 
stars. The simple and quick model allows parameter space ex- 
ploration in detail. Normalization of the piston model shows that 
the normalized 3-D density pattern in the circumstellar envelope 
only depends on two parameters: eccentricity of the star orbit e 
and steady outflow velocity to mean star motion velocity ratio y 
(= V e /u)a p ). When viewed perpendicular to the orbital plane, the 
spiral-like pattern, broken spiral pattern, and incomplete con- 
centric shell pattern can be seen, while when viewed along the 
orbital plane, the alternative concentric half-shell pattern, egg- 
shell pattern, and half-shell half-gap pattern can be seen. Higher 
eccentricity causes more left-right asymmetry and the variation 
of relative radii between the left and right series of concentric 
half shells, while larger y makes a weaker density pattern but 
helps bring out asymmetry from the pattern center. 

One sees a spiral-like pattern in a large observers' inclination 
angle range (0° -79°), while concentric half-shell patterns can be 
seen only in edge-on binary, which indicates that one has more of 
a chance (~ 81%) seeing spiral-like patterns in the sky than see- 
ing concentric shells (~ 19% chance). Fitting Archimedes' spi- 
ral curve to numerical shocked column density patterns demon- 
strates that the spiral center is away from the orbit center by no 
more than 12% of the spiral pitch. The spiral pattern may be- 
come broken when the eccentricity is larger than about 0.4. 
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Appendix A: Formulation of single-piston 
density-pattern 

This section only considers the normalized piston model in 
which the normalized orbital period is 1 and the normalized 
steady outflow velocity is y. Pattern geometry is determined by 
the positions of all density peaks and/or shock regions. Because 
the piston motion is periodic, one only needs to consider gas 
parcels ejected within a single period. Before the occurrence of 
shock, there should be no collision of gas parcels. One may find 
a peak parcel that corresponds to a density maximum. Assuming 
the peak parcel was ejected at an earlier time f p (the peak par- 
cel ejection time), the radial position of the peak parcel (or the 
density peak) at a given time t is 



fpeak 



(t) = l(tp) + [y+Vi(tp)] (f-tp). 



(A.l) 



Here the gas parcel velocity is expressed as two components: 
steady outflow velocity y and an individual velocity component 
V/(f p ) contributed by the star motion. At a later time f c when the 
first collision between gas parcels occurs (the peak parcel is usu- 
ally the parcel that collides with others first), shock region begins 
to form. The time t c is called the shock time, and the parcels that 
collide first are called the "shock parcels''. After that, further gas 
parcel collisions modify the velocity of the shock parcels accord- 
ing to the conservation of momentum. At any given time f > f c , 
assuming the shock parcels were ejected at an earlier time f s (the 
shock parcel ejection time), the radial position of shock parcels 
(or shock region) is 



fshock© = /(fs) + r(f- fs) + ^(fs)(fc - fs) + J V m (t)dt. (A.2) 

The individual velocity component of the shock parcels V m (f) 
in the last term on the right side is different from star motion 
velocity V/(f s ) because collisions modify it according to the con- 
vervation of momentum. V m (t) cannot be analytically expressed, 
but has to be determined by piston model simulation. 

To calculate the radial positions of density peaks or shock 
regions using the above formulas, one needs to know f p ,? s ,? c . 
They can be determined using the definition of density maxi- 
mum. Under the assumption of constant mass-loss rate M and 
outflow velocity y, the overall density distribution p(f, f p ) in the 
circumstellar envelope at a given time t < f c roughly follows an 
inverse square law of radius r(f p ) upon which the density pat- 
terns are superimposed. Therefore, it is helpful to define a pat- 
tern density by multiplying the square of radius by density. It is 
not difficult to deduce the formula of pattern density as 



p'(t, fp) = p{t, fp)r 2 (f p ) = 



M 
4tt 



dV,(t p ) 



(A.3) 
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Fig. 3. Computed maps 
show the azimuth an- 
gle effects on pattern 
geometry. All maps are 
shown for e = 0.8, y = 
10, £ = -90°, that is, 
viewed along the orbit 
plane (edge on), ifr = 
-180° means viewing 
along the major orbit 
axis, while iff = -90° 
means viewing along 
the minor orbit axis. 



When the second inequality in formula ( IA.6I ) becomes equal- 
ity, one can determine the moment when the shock just begins to 
happen, and therefore the solution of the shock parcel ejection 
time f s = f p can be derived by solving formulas ( IA.4IIA.5llA.6l ) 
for this moment. In this case, the shock time can also be derived 
from the second equality of formula ( IA.61 > as 



h + y 



dV,(ts) 
df s 



(A.7) 



With f p , f. 



f c solved from the formulas, one can determine the 
density peak or shock region postions from Eqs. JA.lt and iA.2i , 
and hence the pattern geometry. 

Semi-analytical solutions for the special case of circular or- 
bit have been derived on the basis of the formulas. (The solutions 
are not shown in this paper.) Archimedes' spiral pattern and al- 
ternative concentric shell pattern have been confirmed by the an- 
alytical solutions only in the outer part of the pattern, while, in 
the central region of the patterns, deviation from these standard 
curves is expected. 



By constraining the maximum of the pattern density p'(t, f p ), one 
may derive the following formulas for solving f p : 



d 2 V,(f P ) , v dV,(fp) 



dfn 



d 3 V ; (f P ) ,„ 



df p 3 



< f • 



(f-fp) 



< 2 



dfn 



d 2 V/(f p ) 
dfp 2 



t P <7 



dVKfp) 
df n 



(A.4) 



(A.5) 



(A.6) 
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